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Le but de ce cours est d’apprendre à programmer en fortran 90. Il est prévu 
pour 8 séances d’1h20, accompagné de 8 séances de quatre heures de TD-TP 


sur machine. 


J’ai essayé de le présenter non pas comme une collection exhaustive de com- 
mandes (ce qui permet difficilement de pouvoir programmer avant d’avoir fini 
le cours) mais plutôt comme une progression en trois phases. La première par- 
tie devrait permettre de pouvoir écrire vos premiers programmes après une ou 
deux séances seulement. La deuxiéme introduit certains concepts plus complexes 
(sous-programmes, types dérivées, allocation dynamique). Les troisiéme et qua- 
trième parties mettent l’accent sur les apports réellement nouveaux de fortran 
90 par rapport au vieux fortran 77 (modules, interfaces génériques, surcharge 


d’opérateurs, pointeurs). 


Nombre de commandes présentées ici sont expliquées avec beaucoup plus de 
précision dans divers cours et manuels de référence donnés dans la bibliographie. 
Ces références sont en accés libre et gratuit sur internet. N’hésitez pas a les 


consulter dés que vous en avez besoin. 


Le langage fortran évolue régulièrement, mais la véritable évolution a été faite 
en passant du fortran 77 au fortran 90. Les nouvelles normes fortran 95 (puis 
fortran 2003, fortran 2008) n’apportent que des fonctionnalités nouvelles. Par 


conséquent, tout ce qui est dit ici devrait rester valable avec les normes futures. 


Enfin, je veux insister sur le fait que quel que soit le cours que vous utiliserez, il 
n’y en fait qu’un seul moyen d’apprendre à programmer : c’est de programmer 
soi-même (de la même façon que l’on n’apprend bien une langue étrangère qu’en 


la parlant soi-même). Donc usez et abusez de votre temps en salle machine (le 


temps des séances de TP n’y suffira pas), c’est la garantie du succès. 
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1 Introduction 


1.1 Qu'est-ce qu’un programme infor- 
matique ? 
— tout commence par un algorithme : suite d'opérations 
donnant un résultat final 


exemple : 
approximation de V2 par méthode de Newton 


uo donné 


Un 


1 
calculer Un+1 = I + —, pour n=O a 10. 
Un 


— ces opérations sont écrites dans un fichier, avec un 
langage informatique (pour nous, le fortran 90 (f90)) 
— ce fichier constitue le programme ; 


— c’est une suite d'opérations sur des variables, et 


éventuellement des échanges de données entre uti- 
lisateur et ordinateur ; 

— les variables sont des noms correspondant à de 
la place mémoire dans l’ordinateur, utilisée pour 


stocker les données et les manipuler. 


— exemple : approximation de V2 


program racine 
implicit none 


I--- variables 
integer :: n 
real 


I--- initialisation 


u= 1.0 


Less boucle 
do n=1,10 

u=u /2.0+1.0/u 
end do 


tre affichage 
print*,’approx. de sqrt (2) 


end program racine 


— ce programme n’est pas utilisable tel quel : il faut 
le rendre compréhensible par l’ordinateur, c.-à-d. le 
traduire : le logiciel traducteur est appelé compila- 


teur, la traduction est appelée compilation. 
exemple : 


gfortran -o prog prog.f90 
prog.f90 | : fichier contenant le programme écrit en 


f90 
: traduction (appelé fichier exécutable) 


— le fichier exécutable permet d’exécuter les instruc- 


tion programmées. 


exemple : 


1.2 Le langage F90 : généralités 


— ensemble de règles syntaxiques et grammaticales et 
de mots permettant d'écrire des opérations mathé- 
matiques, de communiquer des informations entre 
l'ordinateur (mémoire), les fichiers (disque dur) et 


l'utilisateur (clavier/écran). 


l'essentiel du cours va porter sur l’apprentissage de 
ce langage : 


— en commençant par une vue rapide et simplifiée 
permettant de pouvoir programmer après une ou 


deux séances ; 


— puis en découvrant ensuite les aspects plus com- 


plexes du langage. 


conventions typographiques utilisées dans ce docu- 
ment : en gras pour les mots-clefs, entre crochets 
[| pour les instructions optionnelles, en italique pour 
les commandes unix, en teletype pour les lignes 


fortran. 


— historique : premier compilateur (langage de haut 
niveau) en 1957, évolutions constantes depuis (voir 
wikipedia) : 

— 1978 : Fortran 77, très répandu aujourd’hui encore dans les 
gros codes de calculs 

— 1990 : Fortran 90, fortran devient un langage moderne (format 
libre, programmation modulaire) 

— fortran 95, 2000 (programmation orientée objet), 2004, ... 


F90 ou C++ pour le calcul scientifique : la majeure 
partie des gros codes de calcul industriel sont écrits 


en fortran ... 


Attention aux effets de mode (“Fortran est dépassé”, 


“vous programmez encore en fortran ?”, “il faut pro- 


grammer en C++”): 
fortran est bien plus simple d’utilisation 
les concepts de C++ ne sont guère utiles pour le calcul scien- 
tifique en général 


les aspects programmation ” orientée objet” du C++ sont main- 
tenant utilisables en fortran (depuis fortran 90 et fortran 2003) 
de nouveaux compilateurs fortran sont développés pour les 
nouvelles architectures (fortran pour les GPU par exemple). 
attention aux commentaires de ceux qui ne connaissent que le 
C (et de ceux qui ne connaissent que le fortran !) 

voir le blog d’un physicien qui connait bien les deux langages : 
http://qd.typepad.com/4/2005/08/c_vs_fortran.html 


2 Premiers pas en F90 


2.1 Exemple de programme 


2.2 structure 


program nom_du_programme 


declaration des variables 


instructions 


[end program nom_du_programme | 


2.3 Les variables 


— variables : noms (chaînes de caractères alphanumériques 
plus quelques autres, cf. section 2.4) permettant de 
manipuler des données en mémoire. 


— opération essentielle : l'affectation (symbole =) 
exemple : variable n (entier) : 


prend la valeur 2 


augmente de 1 


ce n’est pas une égalité mais une opération : la valeur 


de n stockée en mémoire est remplacée par elle-même 
plus 1. 


— le bloc déclaration des variables sert à indiquer à la 


machine quel est le type de chaque variable utilisée 


afin de réserver la place en mémoire nécessaire. 


exemple : 


integer :: n,i | entier n=10 
real :: x,y réel x=1.0e-5 
y=1.0 
complex :: complexe z=cmplx(1.0,2.5e3) 
character :: caractère c=?d? 


logical :: booléen b=.true. 


— déclaration implicite : bien que f90 ne l’exige pas, 
il est prudent de déclarer toutes les variables. Pour 
éviter les oublis, placer l’instruction implicit none 
au début du bloc déclaration. Les oublis sont alors 


détectés à la compilation. 


— les tableaux : tous les types précédents peuvent être 


structurés en tableaux avec l’attribut dimension 


integer, dimension(—2:2) 
real, dimension(1:3,0:4,1:20) 
character, dimension’ i: 3.5476) 


— chaque argument de dimension indique les bornes 
de la dimension correspondante par un intervalle en- 
tier min:max 
— acces aux composantes : 

dim(0)=10 

tabi 2 20) = 2.3e—5 

ton(l.5)= "a" 


2.4 les identificateurs 


caractères du clavier permettant de nommer une va- 
riable ou un autre élément du programme ; 
suite de 1 à 31 caractère alphanumériques parmi : 


lettres sans accents majuscules et minuscules, chiffres, 


caractère _ (souligné) ; 


le 1*’caractère doit être une lettre; 
les majuscules et minuscules ne sont pas différentiées ; 
exemple : 


— les identificateurs suivants sont valides : 


constante_gaz 
pi2 
Rk54 


— les identificateurs suivants ne sont pas valides : 


accentué 

avec espace 
Il_y_a_plus_de_trente_et_un_caracteres 
_souligne_devant 

1_chiffre_devant 

nomalphanumerique 


2.5 Continuation de ligne : caractère & 


une ligne d'instructions comporte au maximum 132 


caractères 


e , . 
si plus de 132 = message d’erreur : 
write(*,*) ’Au clair de la lune, mon ami Pierrot. Prete moi ta plume, pour ecrire un mot. 
Ma chandelle est morte, je n ai plus de feu. Ouvre moi la porte pour 1 amour de dieu’ 

1 


Error : Unterminated character constant beginning at (1) 


solution : couper la ligne en morceaux. Placer & en fin 
de premiere ligne et & au début de la ligne suivante : 


write (*,*x) ’Au clair de la lune, mon ami Pierrot. & 
Prete moi ta plume, pour ecrire un mot. & 
Ma chandelle est morte, je n ai plus de & 


feu. Ouvre moi la porte pour l amour de dieu’ 


2.6 Les commentaires : symbole ! 


tout ce qui est à droite du symbole ! est ignoré par 
le compilateur 

cela permet d'écrire des commentaires dans le pro- 
gramme, c.-a-d. du texte qui en explique le fonction- 
nement 

le ! peut être placé n’importe où dans la ligne (au 
début, au milieu, à la fin), et à n’importe quel endroit 
du programme 


exemple : voir le programme page 30. 


2.7 Les structures de contrôle 


Ce sont des instructions qui permettent de faire des 


itérations ou des tests sur des variables. 


2.7.1 Le test if 


— exemple : calcul de la valeur y = { 7'97 si 770 
sinon 


if (x/=0.0) then 
y = xe lop (x) 
else 
y = 0.0 
end if 


— syntaxe générale : 


[nom:| if (expression logique) then 
instructions 

|elseif (expression logique) then 
instructions | 


| else 
instructions | 
end if [nom] 


où nom permet de donner un nom à la structure pour 


une meilleure lisibilité du programme (par exemple 


en cas de if imbriqués) 


— expression logique est une combinaison d’un ou 


plusieurs tests sur des variables, ainsi que des trois 


opérateurs logiques 


égal 

différent 
supérieur strict 
inférieur strict 
supérieur ou égal 
inférieur ou égal 
et 


ou 


— les variables de type logical peuvent aussi être uti- 
lisées 
— le else if et le else sont optionnels 


O 31 27<0 on 572 
— autre exemple : calcul de y = 4 Vz si xe[0,1] 
1 


si x>l 


rac:if (x<0.0.or.x>2.0) then 
y = 0.0 
elseif (x<=1.0) then 


y=sqrt (x) 
else 

y=1.0 
end if rac 


2.7.2 La sélection select case 


— pour choisir une instruction à exécuter en fonction 


de la valeur d’une expression scalaire de type entier, 


caractère ou logique 


— exemple : 
character(len—30) :: pays 


langue : select case(pays) 
case(’france’,’quebec’,’suisse 
printx, ’? bonjour’ 
case( ’royaume—uni’,’usa’) 
printx,’ hello’ 
case default 
printx, langue pas disponible 
end select langue 


1 


,’ belgique’) 


1 


— l’expression testée doit vérifier au plus un case 


— le case default n’est pris en compte que si l’expres- 


Sion ne vérifie aucun case 


2.7.3 Itérations : structure do 


— forme énumérative (nombre d’itérations fixé) 


ni 
exemple : calcul de y = > _ 
1 


i=1 


integer :: i,n 
real :: x,y 


y —0.0 
do i=Il1,n 


y=y+x*xi/i 
end do 


syntaxe : 


[nom:|do variable=debut ,fin |, pas |] 
instructions 
end do [non] 


autre exemple : calcul de s = 1 +3+5+7+9 
integer 


s=0 

do i=1,9,2 
s=s+i 

end do 


— forme infinie : on ne sort qu'avec l’instruction exit 


[nom:]| do 
instructions 
if (condition) then 
exit [nom] 
end if 
instructions 
end do [nom] 


exemple : 


do 
read * ,nombre 
if (nombre==0) then 
exit 
end if 


somme — somme + nombre 


end do 


— forme conditionnelle : arrêt si la condition spécifiée 


en début de boucle est satisfaite 


[nom:] do while (condition) 
instructions 
end do [nom| 


exemple : 


reponse=’oul’ 

lecture:do while ( reponse’ oui’) 
printx,’continuer ?’ 
read * ,reponse 

end do lecture 


2.8 Les entrées-sorties écran/clavier 


— lecture au clavier de la valeur d’une variable, quel 
que soit son type 


readx,variable 


— écriture à l’écran d’un message et/ou de la valeur 
d’une variable 
print x,’ bonjour’ 
print x, x 


print,’ il est’ ,h,’ heures’ 


— exemple : 


integer :: n 
character(len=15) :: pays 


printx,’ quel age avez-vous ?’ 
read x ,n 

print*x, vous avez 
printx,’de quel pays etes vous ?’ 
read x ,pays 

prints,’ vous venez de ’ ,pays 


1 


2.9 Compilation 


pour faire fonctionner un programme écrit en f90, il 
faut le traduire en un langage compréhensible par 


l’ordinateur : c’est la compilation. 


cette compilation est effectuée par un compilateur. 
Dans ce cours, nous utiliserons le compilateur libre 
gfortran. 

la compilation crée un fichier dit exécutable. 
exemple : le programme s’appelle calcul, il est écrit 
dans le fichier toto.f90. On souhaite compiler ce pro- 
gramme pour créer un exécutable nommé ici com- 


pute. On tape donc la commande 


gfortran -o compute toto.f90 


plus généralement, la commande de compilation est : 


gfortran -o executable fichier_fortran 


ou executable est le nom que l’on souhaite donner au 


fichier exécutable et fichier_fortran est le nom du fi- 


chier qui contient le programme f90 (il est conseillé 


de donner au fichier le même nom que le programme, 


en rajoutant l'extension .f90). 


— différentes options existent, cf. section 3. 


— messages d’erreurs : 


— le compilateur détecte certaines erreurs (syntaxe, 


conformité de types, de dimensions) 


il affiche un message donnant le numéro de la ligne 
contenant l'erreur, avec une brève explication en 
anglais 


exemple : (programme contenu dans le fichier erreur_exemple.f90) 


erreur_exemple.f90:7.3: 


y = x*log(x) 
1 
Error: Symbol ’y’ at (1) has no IMPLICIT type 


une erreur en provoque souvent d’autres (exemple : 


une variable non déclarée va générer des erreurs a 


chaque ligne ot elle est utilisée) 


la liste des erreurs est donc souvent tres longue, et 
il faut alors commencer par corriger les premieres 


de la liste, puis recompiler. 


— les messages d'erreurs ne sont pas toujours tres 


clairs ... 


exemple : dans un programme contenu dans le fi- 


chier erreur_declaration.f90, on fait une erreur de 
frappe dans la déclaration suivante : 
ineger :: n 
au lieu de : 
integer :: n 
ce qui provoque le message suivant : 
erreur_declaration.f90:5.2: 
ineger :: n 


1 


Error: Unclassifiable statement at (1) 


2.10 Exécution 


si le fichier contenant l’exécutable s’appelle [prog | 


alors taper simplement 


erreurs à l’exécution 


— toutes les erreurs ne sont pas détectées a la com- 
pilation (certaines options peuvent cependant y 
aider) 


— divisions par 0, erreurs d’affectation de tableau 
— erreurs d’algorithme (+ au lieu de —) 


— il faudra apprendre a trouver ces erreurs en fai- 


sant des tests, et en utilisant un debugger (cf. 
sec. 3.5.5). 


2.11 Ecrire un programme proprement 


— écrire des commentaires (but et fonctionnement du 
programme, signification des variables, décrire chaque 
bloc d'instruction) 
indenter correctement (en particulier les structures 
de contrôle) 
nommer les boucles imbriquées 
déclarer les variables dans l’ordre croissant des struc- 
tures (logiques, entiers, caractères, réels, puis tableaux) 
mieux vaut utiliser des noms de variables longs mais 


clairs que des noms courts et obscurs 


il est fortement conseillé d'utiliser l’éditeur de texte 
Emacs pour programmer. Outre sa puissance d’édition 
en général (au moyen de raccourcis clavier qui per- 
mettent de taper du texte très rapidement), il possède 
des fonctions adaptées à la programmation : inden- 
tation automatique, complétion des noms de pro- 
grammes et de boucles, coloration des mots clefs, 
etc. Voir le guide Emacs distribué en cours [5]. 


2.12 Exemple complet 


£ 


calcul de la solution de x — e7” = 0 par la méthode de 


Newton. 
algorithme : 
zo et tol donnes 


ee) 


caleuler Tagi = in pue 
n41 — T 
tant que paa > tol 
n 


program zero 
implicit none 


l--- variables 
real :: x,y,diff ,tol 


'--- lecture des données 
print», entrer le x initial’ 
read *,x 


printx,’tolérance ? 


read x,tol 


!'--- méthode de newton 
newton:do 


'-- ~tération 


y=x—(x—exp(—x))/(1.0+exp(—x)) 


!-- différence entre deux 1itérés 


if (x/=0.0) then 
diff=abs ((x—-y)/x) 
else 
diff=abs(x-—y) 
if 


l-- mise à jour 
x=y 


!-- test de convergence 
if (diff<tol) then 


exit 
end if 


end do newton 
print,’ solution approchée=’ ,x 


end program zero 


3 Un peu plus loin ... 


3.1 Variables “constantes” 


Pour déclarer une variable “constante” (c.-à-d. qui a 


une valeur non modifiable) : 

— attribut parameter 
la variable est non modifiable, sinon erreur détectée 
à la compilation 
utile pour déclarer et protéger des constantes phy- 
siques 
exemple : 
real, parameter :: pi=3.14159265, N=6.02252e23 


dans la déclaration, les opérations usuelles sont uti- 
lisables, mais pas l’appel de fonction : 

real, parameter :: quart=1.0/4.0 !--- valide 

real, parameter :: pi=acos(—1.0) /--- non valide 


3.2 Sous-programmes 


3.2.1 Objectif 


— bloc d'instructions utilisé plusieurs fois = l’écrire 


dans un sous-programme une fois pour toutes, et 


l’appeler quand nécessaire (aussi appelé procédure) 


— deux types de sous-programme : les subroutines et 


les fonctions 


— exemple 1 : subroutine 


program exproc 
implicit none 


real :: moyenne, maximum 
real, dimension(100) :: tab 


call random_number(tab) 
call sp( tab, moyenne, maximum ) 
print *,moyenne , maximum 


contains 
subroutine sp( t, moy, max ) 
implicit none 


real, dimension(100), intent(in) 
real, intent(out) :: moy, max 
integer :: i 


max = t( 1); moy = t(1) 
do i=2,100 
if (t(i) > max) then 
max = t(i) 
end if 
moy = moy + t(i) 
end do 
moy = moy/100 


end subroutine sp 


end program exproc 


— exemple 2 : fonction 
program exfct 
implicit none 


real :: maximum 
real, dimension(100) 


call random_number(tab) 


maximum — maxi(tab) 
print *,maximum 


contains 
function maxi(t) 


implicit none 


real, dimension(100), intent(in) 
integer :: i 
real :: maxi 1- La fonction est déclarée 


maxi = t(1) 
do i=2,100 
if ( t(i) > maxi ) then 
maxi = ti) 
endif 
end do 
end function maxi 


end program exfct 


3.2.2 Fonctionnement 


— place dans le programme : 
avant la dernière ligne end program nom_programme 
et après le mot-clef contains 
program nom_programme 
declarations 
instructions 


contains 
subroutine nom subroutine 


end subroutine nom subroutine 


function nom _ fonction 


end function nom_ fonction 
end program nom_programme 


on parle alors de sous-programme interne (par oppo- 


sition aux sous-programmes externes, cf. section 4.3) 
portée des variables : toutes les variables du pro- 
gramme principal sont connues de ses sous-programmes 
internes. Pas besoin de redéclaration, sauf pour les 
arguments (voir plus loin). 


les variables utilisées seulement par un sous-programme 
(variables locales) : déclarées à l’intérieur de celui-ci, 
comme dans un programme principal. Connues de 
lui seul. 


— appel d’un sous-programme : 


— subroutine : appelée par l'instruction call 


call nom _subroutine (arguments) 


— fonction : appelée comme une fonction mathé- 


matique : 


y=nom_fonction (arguments ) 


3.2.3 Arguments d’une subroutine 


subroutine : sous-programme ayant ou non des argu- 


ments d’entrée et de sortie : 


— variables à déclarer à l’intérieur de la subroutine 


— il est conseillé de préciser dans ces déclarations les 
attributs suivants 
intent(in) pour les arguments d’entrée 
intent(out) pour les arguments de sortie 
intent(inout) pour les arguments mixtes. Ainsi le 


compilateur détecte une mauvaise utilisation des ar- 


guments (entrée modifiée, sortie non donnée, etc.) 


3.2.4 Arguments d’une fonction - déclaration 


du résultat 


fonction : comme une subroutine, mais 
— arguments d'entrée obligatoires 
— renvoie nécessairement un résultat, stocké dans une 


variable ayant le même nom que la fonction 


— le type du résultat peut être donné soit avant le nom 
de la fonction (sauf si résultat tableau) : 


real function maxi(t) 


end function maxi 


soit dans une déclaration : 
function maxi(t) 
veal :: maxi 
ed: function maxi 


— le résultat peut avoir un nom différent du nom de la 


fonction, en utilisant la clause result 


function maxi(t) result(y) 


real :: y 


end function maxi 


dans ce cas, pas de type avant le mot clef function. 
— la variable résultat ne doit pas avoir d’attribut intent 


3.2.5 Les tableaux en arguments 


— problème : dans l’exemple de la page 34, la subrou- 
tine sp n'accepte que des tableaux de taille 100. 


Comment faire pour qu’elle fonctionne avec des ta- 


bleaux de tailles quelconques ? 


— le plus simple : donner en entrée le tableau et ses 


dimensions (”assumed size array” en anglais) : 


program exproc 
real, dimension (100) 


contains 
subroutine sp( t, n, moy, max ) 


integer , intent(in) :: n 

real, dimension(l:n), intent(in) 
real, intent(out) :: moy, max 
integer :: i 


max = t( 1); moy = t(1) 
do i=2,n 


end do 
moy = moy/100 


end subroutine sp 


end program exproc 


— autre exemple : 


program toto 


real, dimension(1:3 ,1:2) 


call sousprog(A,3,2) 


contains 
subroutine sousprog(A,m,n) 
integer, intent(in) :: m,n 
real, dimension(1:m,1:n) 


end subroutine sousprog 
end program toto 


— on peut aussi ne passer que le tableau en argument 
et récupérer ensuite ses dimensions (tableau à profil 
implicite, ” assumed shape array” en anglais). Atten- 


tion alors a la déclaration : 


program exproc 


real, dimension (100) 


contains 
subroutine sp( t, moy, max ) 


real, dimension(:), intent (in) 
real, intent(out) :: moy, max 
integer :: i,n 

n=size(t) 
max = t(l 
do i=2,n 


); moy = t(1) 
end do 
moy = moy/100 

end subroutine sp 


end program exproc 


— autre exemple : 
program toto 


real, dimension(1:3,1:2) :: A,B 
real, dimension(1:3) :: u,v 


call sousprog(A,B) 
y=fct (u) 


contains 

subroutine sousprog(A,B) 
real, dimension (:, 
real, dimension (:, 
integer :: m,n 
m=size(A,l) ; n=size(A,2) 
B=A * *2 

end subroutine sousprog 


intent (in) 
intent (out ) 


1 


:) 
:) 


? 


function fct(u) result(v) 
real, dimension(:), intent (in) 
real, dimension(size(u)) :: v 


end function fct 
end program toto 


— ATTENTION : le profil implicite s’utilise seulement 
pour les tableaux arguments. Pour des tableaux lo- 
caux ou un résultat de fonction dont les tailles dé- 
pendent de celles des arguments, remarquer l’utili- 


sation de l’instruction size. 


3.2.6 Les chaînes en arguments 


même principe que pour les tableaux : 


— on peut passer en argument la chaîne et sa longueur 


— on peut aussi utiliser un profil implicite, pour les 


arguments uniquement 


— exemple : 


program chaine 


character(len—20) 


printx,ch,f2(ch,20) 
print+,f1(ch) 


contains 


function f1(ch) 
implicit none 
character(len=x), intent(in) :: 
character(len=len(ch))) :: f1 


end function eff_blanc 
function f2(ch,n) 
implicit none 
integer, intent(in) :: n 
character(len=n), intent (in) 
integer :: £2 


end function nb_blancs 


end program chaine 


3.3 Tableaux 


3.3.1 Tableaux statiques 


déjà vu : les dimensions sont fixées à la déclaration : 
real, dimension(1:10,1:4) :: A 


3.3.2 Allocation dynamique 


— il est fréquent que les tailles des tableaux utilisés 


dans un programme dépendent de données et ne 
soient pas connues a la compilation 

moyen simple et efficace : attribut allocatable 
(tableau allouable) : 


real, dimension(: ,:), allocatable 


— le nombre de dimensions est donné (ici 2), mais pas 
les dimensions elles-memes 
— pour manipuler ces tableaux, on commence par les 


allouer en mémoire (c.-à-d. la place nécessaire est 


réservée), en donnant leurs dimensions. Exemple : 


integer 


read*x,n,m 
allocate(A(1:n,1:m)) 


— pour changer la taille du tableau : le désallouer, puis 
le réallouer : 


deallocate(A) 
allocate(A(1:3,0:4)) 


— évidemment, on ne peut utiliser un tableau non en- 
core alloué. 


ne pas oublier de libérer la place mémoire avec deallocate 
quand le tableau n’est plus utilisé (en particulier 
pour les tableaux locaux dans des sous-programmes, 


sinon, risque de dépassement de mémoire). 


3.3.3 Allocation dynamique dans un 


sous-programme : tableaux automatiques 


autre moyen simple de créer un tableau dont la taille 
n’est pas connue à la compilation 


utilisable seulement au sein d’un sous-programme 
principe : les dimensions du tableau sont des va- 
riables passées en arguments du sous-programme (” au- 
tomatic array” en anglais). Exemple : 


subroutine autom(m,n) 


integer :: m,n 
real, dimension(m,n) 


avantage : plus rapide en exécution que la création 
d’un tableau allouable. A utiliser pour les sous-program- 
mes appelés très fréquemment. Les instructions for- 
tran sont aussi plus simples à écrire. 

inconvénient : aucun moyen de vérifier s’il y a assez 

de place pour créer le tableau, contrairement aux 


tableaux allouables. Le message d’erreur correspon- 


dant est généralement très obscur. 
pour les tableaux de grandes tailles utilisés peu fréquem- 
ment, il est donc plutôt conseillé d’utiliser les ta- 


bleaux allouables. 


3.3.4 Terminologie 


rang : nombre de dimensions (le rang de A est 2) 


taille d’une dimension : nombre d’éléments dans cette 
dimension (les tailles des dimensions de A sont 3 et 
5) 

forme : suite ordonnée des tailles des dimensions (la 
forme de A est 3,5) 

une taille peut être nulle 


3.3.5 Opération conformes entre tableaux 


— but : faire des opérations sur des tableaux sans écrire 


les opérations élément par élément 


deux tableaux sont conformes s’ils ont même forme : 
A(1:10,1:4)et B(3:12,0:3) sont conformes 
A(1:10,1:4) et B(3:12,1:3) ne sont pas conformes 
A(1:1,1:1) et B(1:1) ne sont pas conformes 


une expression est conforme si elle fait intervenir des 
tableaux conformes 

un scalaire est considéré comme conforme à tout ta- 
bleau 

la plupart des opérations scalaires sont applicables 
aux tableaux : opération appliquée élément par élément, 
résultat de même forme que les tableaux opérandes 


— exemples : 


program conforme 


implicit none 


integer :: i,j 
real, dimension(3,3) :: A,B,C 


call random_number (A) 
call random_number (B) 


C=A+B+2.0 !-- op. équivalente à 
print x,C 
do i=1,3 
do j=1,3 
C(i,j)=A(i,j)+B(i,j)+2.0 
end do 
end do 
print x,C 


C=A*B |== pas prod, matricrel 
print x ,C 

C=sqrt(B)*xC 

print x,C 

C=B+cos(Bx3.1416) 

print x ,C 

A=1—sin(log(B+A)) 

printx,A 


end program conforme 


3.3.6 Créer des tableaux “a la main” 


— avec une liste de variables d’un méme type 
vari, var2,..., varn, on peut créer un vecteur (ta- 


bleau de rang 1) avec le constructeur (/.../) : 


v=(/vari,var2,...,varn/) 


— la liste peut-elle même être créée avec une boucle 
implicite 
hr La) !--équivaut à 
=(/1,4,9,16/) 


— on peut imbriquer plusieurs boucles implicites 


=(/((itj J=] l -équivaut à 
=(/2,3,3,4,4,5/) 


— comme dans une boucle do on peut donner un pas : 


— pour transformer un vecteur en tableau à plusieurs 
dimensions : 
— fonction reshape(vecteur,forme_du_tableau) 


real, dimension(2,3) 
A=reshape(v,(/2,3/)) 


donne la matrice A = (42 f ) (remplie colonne par 


colonne) 


3.3.7 Les sections de tableaux 


on accède à une section de tableau en modifiant conve- 


nablement les indices de début et de fin des dimen- 


A 
Lignes de 2 à 6 
idem 
ieme ligne 
jeme colonne 
élément communs aurt 
lignes 2 à 6 et 


colonnes 1 à 4 


3.3.8 Stockage des tableaux dans la mémoire 
et ordre des boucles 


— il est très important lorsqu'on manipule de gros ta- 


bleaux de savoir comment ils sont stockés en mémoire 


— exemple : y-a-t’il une différence importante entre les 


deux blocs suivants ? 


=> oui : le premier bloc est beaucoup plus coûteux 


en temps! 


— quel que soit le nombre de dimensions d’un tableau. 
il est stocké en mémoire sous forme uni-dimensionnelle. 


Les éléments sont stockés ainsi : 


an aan eo] [en [an ETC ETC 


— exemple : les éléments du tableau suivant 


integer, dimension(2,4) :: a 
sont rangés ainsi : 
a(1,1) a(2,1) a(1,2) a(2,2) a(1,3) a(2,3) a(1,4) a(2,4) 


— conclusion : quand on parcourt les éléments d’un ta- 


bleau, il faut que la boucle sur la première dimension 


soit la boucle la plus interne (ceci est lié à la notion 


de cache dans la mémoire) 


— c’est l’inverse de ce qui se fait en C, et de ce qui se 
fait en général quand on écrit les algorithmes 
— voir en TP pour une comparaison en temps calcul 


des deux stratégies. 


3.4 Entrées-sorties (E/S) 


On a vu dans la première partie comment lire au clavier 
et écrire à l’écran. Ici, on va voir comment lire et écrire 


dans un fichier. 


3.4.1 E/S en format libre, dans un fichier texte 


à accès séquentiel 


C’est le plus simple, équivalent de ce qu’on a vu avec 


l’écran/clavier 
— ouverture du fichier toto : 


open(unit=10, file=’ toto’) 


— cette instruction ouvre le fichier à la première 
ligne (s’il n’existe pas, il est créé). 
— | expression unit=10 indique que le fichier toto 


est connecté au descripteur numéro 10. 


— tous les nombres peuvent étre utilisés comme 


descripteurs sauf 5 (clavier) et 6 (écran). 


— si le mot clef file n’est pas donné, le fichier 


fort.n| est créé, ou n est le descripteur donné 


après le mot clef unit. 


— écriture de données dans toto (si associé au descrip- 
teur 10) : 


write(10,*) donneel, donnee2, donnee3 


— une instruction write écrit toutes les données sur 
la méme ligne, le write suivant passera a la ligne 


Suivante. 


— les données peuvent étre de n’importe quel type 


(constantes ou variables) 


— le print*, pour écrire à l'écran est équivalent a 


write(6,*) et à write(*,*) 


— lecture de données dans toto (si associé au descrip- 
teur 10) : 


read(10,*) vari, var2, var3 


— la premiere ligne de toto est lue, le read suivant 


passera a la ligne suivante. 


les données sont supposées étre séparées par des 
blancs. 


elles sont stockées respectivement dans vari, var2, 


.., Varn. 


si plus de données que de variables : seules les 


premières données sont stockées (cas inverse : seules 


les premières variables sont affectées) 


le readx, pour lire au clavier est équivalent à 
read(5,x) et à read(x,*). 


— fermeture du fichier toto : 


close(10) 


— attention, si on réouvre le fichier avec open, on se 


retrouve à la première ligne. 


— écrasement de fichier : un write après l’ouverture 


d’un fichier existant détruit le contenu du fichier. 


— exemple 1 : 
fichier titi.txt : 


bonjour 4 heures 
temperature See 


programme : 


character(len=14) :: c1,c2,c3 
integer :: n 
real :: T 


open(unit=10,file=’titi.txt’) 

read(10,*) c1,n,c2 !-- tecture lere ligne 
read(10,*«) c3,T pes 2eme ligne 
close (10) 

open( unit =20, file=’toto.txt’) 


!-- ecriture lere ligne 
write(20,*) c1i,’il est ’, n,c2 


!/-- ecriture 2eme ligne 
write(20,*) ’la ’,c3,’est de ’,T,’degres’ 


close (20) 


apres exécution, fichier toto.txt : 


bonjour il est 4 heures 
la temperature est de 37.20000 degres 


— exemple 2 : dans le fichier info.txt, liste de 1000 per- 


sonnes avec un nom, age et taille sur chaque ligne. 


Programme de lecture : 


integer :: i 

character (len=30), dimension(1:1000) 
integer , dimension(1:1000) :: age 
real, dimension(1:1000) :: taille 


open( unit =30, file=’info.txt’) 
do i=1,1000 
read(30,*) nom(i), age(i), taille(i) 
end do 
close (30) 


3.4.2 E/S dans un fichier binaire à accès 
séquentiel 


— les fichiers de textes sont des fichiers de caractères, ce 
qui peut prendre beaucoup de place s’il y a beaucoup 
de données. 

— dans ce cas, il vaut mieux stocker les données (en 
particulier numériques) sous forme binaire. Le fichier 
n’est pas lisible avec un éditeur de texte, mais il est 
beaucoup plus petit. 

— ouverture : avec open comme pour un fichier texte 


en précisant le paramètre form=’unformatted’ : 


— lecture/écriture : open et write sans le paramètre * : 


read(15) vari, var2 
write(15) donnee 


— exemple : tableau de réels aléatoires, stocké dans un 


fichier texte et dans un fichier binaire. Comparer les 


tailles des deux fichiers. 
real, dimension(1:1000) :: tableau 
call random_number (tableau) 


open(unit=1,file=’texte.txt’) 
write(1,*) tableau 
close (1) 


open( unit =2,file=’ binaire.b’ ,form=’ unformatted ’ ) 
write(2) tableau 
close (2) 


3.4.3 E/S formattées 


— inconvénient de E/S en format libre : 


— fichier pas propre : réels écrits sous forme décimale 


si > 10772, sous forme scientifique sinon (cf le fi- 


chier texte.txt de la section 3.4.2). 


— lecture peu sure : risque de mélange des variables 


— format : on indique précisément le nombre et le type 
des variables à lire et écrire 
— paramètre fmt de read et write : 


read(descripteur ,fmt=’(chaine)’) vari ,... 
write(descripteur ,fmt=’(chaine)’) vari ,... 


où chaine indique le type des données dans l’ordre 
où elles sont lues ou écrites 
— formats de base (plus de détails sur [3]) : 

— fmt=’(nIm)’ : n entiers de m chiffres 

— fmt=’(nFm.d)’: n réels sous forme décimale, formés 
de m caractères (dont le .), dont d après la virgule 

— fmt=’(nEm.d)’:n réels sous forme scientifique, formé 
de m caractères, dont d après la virgule (dont le .) 

— fmt=’(nAm)’ : n chaînes de m caractères 

Dans tous les cas, c’est le caractère espace qui sépare 

les n champs. 


— exemple : 


character(len—10) :: prenom, nom 
integer :: age 


prenom=’JEAN ’ 
nom= DURAND’ 
age —20 


write(*x,fmt—’(2A10,113)°) prenom, nom, age 


va écrire : 


— si on doit utiliser plusieurs fois un même format, 


on peut stocker la chaîne associée au paramètre fmt 


dans une variable. Exemple : 


formate—’(2A10,113) 
write(*,fmt-formate) prenom, nom, age 


3.4.4 Contrôle des erreurs en E/S 


avec les instructions précédentes, le programme s’arrêt 


à la moindre erreur (lecture dans fichier inexistant, 
lecture dépassant la taille du fichier, écriture texte 
dans un fichier binaire, etc.) 

— des paramètres des commandes open, read et write 


permettent d'éviter l’arrêt et de déterminer la source 
de l’erreur. 


— exemple : le paramètre tostat : 


integer :: erreur 


read(15,*,iostat=erreur) !- fich. texte 
read(23,iostat=erreur) !- fich. binaire 


— si lecture correcte => erreur vaut 0 


— sinon, si erreur < 0 alors fin de fichier ren- 


contrée, sinon autre erreur. 


— exemple: 


erreur=0 

do while (erreur==0) 
read(15,*x,iostat—erreur) var 
print*+,’/variable=”,var 


end do 


if (erreur <0) then 


1 


printx,’fin du fichier 
else 
print», erreur de lecture 


end if 


1 


— voir les instructions open et write pour la signi- 


fication de iostat. 


3.5 Trouver des erreurs dans un pro- 


gramme 
3.5.1 Les éviter et les détecter 


. avant de programmer un algorithme, s’assurer qu’il est 


correct 


. programmer proprement (voir sec. 2.11), utiliser au 


maximum les signaux d’erreurs renvoyés par les com- 


mandes (paramètre iostat, etc.), afficher des messages 
d'erreur en cas d'opération interdite (x < 0, etc.) 


. lors de la phase de développement, utiliser des options 


de compilation pour détecter les erreurs à l’exécution 


comme | -fcheck=all | qui détecte en particulier les dé- 
passements de tableau) ou| -ffpe-trap=invalid,zero | qui 


détecte les opérations invalides 


. effectuer des tests judicieux pour tester le programme 


. en désespoir de cause, rajouter des affichages de va- 
riables dans le programme, et utiliser l’instruction stop 


pour que le programme s’arrête à l’endroit voulu. 


La tâche 5 peut être simplifiée en utilisant un “debugger” 
(voir sec. 3.5.5). 


3.5.2 Erreurs détectées à la compilation 


— le compilateur renvoie un message plus ou moins ex- 


plicite (voir sec. 2.9) 


— erreur de syntaxe : exemple la compilation sur natty 


de 


program erreur_syntaxe 
implicit none 


integer 
real 


r=10 2 
end program erreur_syntaxe 


renvoie les messages suivants : 


erreur_syntaxe.f90:4.11: 


integer : x 
{ 
Error: Invalid character in name at (1) 


erreur_syntaxe.f90:7.2: 


r=10"2 
1 


Error: Unclassifiable statement at (1) 


— oubli de déclaration : voir l’exemple précédent 


— erreur de dimension : la compilation sur natty de 


program erreur _dimension 
implicit none 


real, dimension(1:4) 


a(4,4)=1.0 
end program erreur_dimension 


renvoie les messages suivants 


erreur_dimension.f90:6.3: 


a(4,4)=1.0 
1 


Error: Rank mismatch in array reference at (1) (2/1) 


— malheureusement, il existe beaucoup d’autres erreurs 
possibles ... 


3.5.3 Erreurs détectées à l’exécution 


Il y a évidemment beaucoup d'erreurs possibles. En 


voici quelques unes avec le programme suivant 


program erreur _execution 
implicit none 
real :: x 


x=1/3 
print x*,x 


x=—1.0 
print *,sqrt(x) 


read x« ,x 
end program erreur_execution 


En supposant que l'utilisateur tape la lettre a au cla- 


vier, on aura alors les affichages suivants a l’écran : 


0.0000000 
NaN 
a 
At line 11 of file erreur_execution.f90 (unit = 5, file = ’stdin’) 


Fortran runtime error: Bad real number in item 1 of list input 


le premier print renvoie 0, car la division est entière 


(oubli de .0 sur le numérateur et/ou dénominateur). 


Mais il n’y a pas d’erreur d’exécution à proprement 


parler, puisque le programme continue a tourner) 


le deuxiéme renvoie NaN (not a number), car la racine 
d’un nombre négatif n’est pas définie (la non plus, 
pas vraiment de message d’erreur : le programme 


continue à tourner) 


le dernier affichage est un message d’erreur : le pro- 
gramme s’arréte (on dit qu’il “plante” ) à cause de la 


lecture d’un caractère au lieu du réel attendu. 


3.5.4 Erreurs d’algorithme 


Ce sont des erreurs souvent difficiles à trouver. C’est 
pourquoi il vaut mieux bien écrire son algorithme avant 


de le programmer 


3.5.5 Les “debuggers” 


La tâche 4 décrite à la section 3.5.1 peut être simplifiée 


avec les logiciels dits “debuggers” : 


— sans eux, il faut recompiler et exécuter le programme 
après chaque modification, ce qui peut être long quand 


le programme est gros ; 


un debugger permet, sans modifier le programme, 
d'afficher la valeur des variables au cours de l’exécution, 
et de stopper temporairement l’exécution à l’endroit 
désiré. 
l’utilisation d’un debugger nécessite un certain inves- 
tissement, mais cet effort est largement récompensé 
par la suite. 

Sur la machine natty, le debugger est gdb. Il peut être 


utilisé avec l’interface graphique ddd (voir le manuel [2] 


avec un exemple simple et instructif dans la section 


” sample session” ). 


Important : pour utiliser un debugger, il faut compiler 
le programme avec l’option 1-8 | 


4 Programmation modulaire 


4.1 Introduction 


4.1.1 Qu’est-ce-que c’est ? 


— programmer de facon modulaire consiste a découper 
le plus possible un programme en morceaux algorith- 
miquement indépendants (appelés sous-programmes) 
exemple : calculer t = trace(A~!B). Ce problème 


contient trois calculs qui peuvent être programmés 
indépendamment : 


A + trace(A) 
(A,B)+ AB 
Am A? 


on essaiera donc de rendre ces morceaux les plus 


indépendants possible au niveau informatique. 


4.1.2 Intérêt 


permet de réutiliser plusieurs fois un morceau de pro- 
gramme facilement 

le programme est plus lisible 

le risque d’erreurs de programmation est moins grand : 


chaque bloc est écrit, compilé, et testé séparément. 


exemple : création d’un ensemble de fonctions pour 


faire des calculs matriciels (produit matrice-vecteur, 


triangularisation, solution d’un système linéaire, dé- 


terminant, trace, etc.) 


on a déjà vu un moyen de créer des sous-programmes 
en utilisant les fonctions et subroutines. On va aussi 


voir les modules. 


4.2 Les procédures internes 


déjà vu à la section 3.2, facile à faire 
problème : pas emboîtable. Une procédure interne 
ne peut pas contenir de sous-programme. Or une 


procédure interne peut nécessiter l’utilisation d’une 


autre procédure interne (exemple : la procédure de 


résolution d’un système linéaire par méthode de Gauss 


nécessite une procédure interne de résolution d’un 
système triangulaire) 

problème de lisibilité : le programme principal peut 
être très gros car il contient touts les sous-program- 


mes 


cette méthode ne répond donc pas parfaitement au 


problème. 


4.3 Les procédures externes 


4.3.1 Pourquoi “externes ?” 


Ces procédures sont écrites à l’extérieur du programme 

principal : 

— soit après l'instruction end program qui termine le 
programme principal 


— soit dans un autre fichier 


4.3.2 Utilisation : probleme d’interface 


Elles peuvent s’utiliser exactement comme une procédure 
interne (même syntaxe), mais il existe un problème 
d'interface : 

— programme principal et procédures sont compilés sé- 


parément 


— il n’y a pas de contrôle de cohérence des arguments 
entre l’appel et la procédure appelée : les erreurs 
éventuelles ne sont donc pas détectées à la compi- 
lation. Elles seront “visibles” à l’exécution, mais se- 
ront peut être difficiles à localiser. Exemple : types 
des variables différents entre l’appel et la procédure 


appelée. 


on dit que l’interface entre programme principal et 
procédure appelée est implicite (contrairement aux 


procédures internes ou l'interface est dite explicite). 


conséquence : il est impossible de passer des tableaux 
ou des chaînes à profils implicites dans une procédure 


externe avec interface implicite. 


=> pour rendre une interface explicite, le plus simple 


est d’enfermer la procédure externe dans un module 


(il existe d’autres solutions moins simples avec les 


bloc interface, voir [3]). 


4.4 Les modules 


— utilisés pour que l'interface entre programme princi- 


pal et procédure externe soit explicite 


— on enferme la procédure externe dans un module de 


la facon suivante : 


module nom_du_module 
implicit none 
contains 


subroutine nom_subroutine 


end subroutine nom_subroutine 


end module nom_du_module 


— le module est alors placé dans un fichier (différent 


de celui du programme principal) qui a pour nom 
nom_du_module.f90 

— dans le programme utilisant la subroutine, il suffit 
de rajouter l'instruction use nom_du_module en tout 


début de programme, avant l’instruction implicit none. 


— exemple : erreur d’argument dans l’appel d’une pro- 
cédure externe 


1. interface implicite (sans module) 
fichier test1.f90 


program testi 
implicit none 
call somme(’oui’,’non’) 


end program testi 


fichier somme.f90 


subroutine somme(x,y) 
implicit none 


real, intent(in) :: x,y 
real :: s 


S=x+y 
print x, somme—’ ,s 


end subroutine somme 


en entrée : 2 chaines au lieu de 2 réels => erreur pas 


détectée a la compilation. A l’exécution, on a un 
résultat aberrant. 


2. interface explicite avec module 
fichier test2.f90 


program test2 
use mod_somme 
implicit none 


"y non’) 


call somme (’ oui 


end program test2 


fichier mod_somme.f90 


module mod_somme 
implicit none 
contains 
subroutine somme(x,y) 
implicit none 


real, intent(in) 
real :: s 


S=x-+y 
print x, ’somme—’,s 


end subroutine somme 


end module mod_somme 


l'erreur est alors détectée à la compilation : 


test2.f90:7.13: 


call somme(’oui’,’non’) 
1 
Error: Type mismatch in argument ’x’ at (1); passed CHARACTER(1) to REAL(4) 


— autres utilisations des modules 


— permettent de définir des constantes, des variables, 


des types dérivés que l’on peut utiliser dans différents 
sous-programmes sans les passer en argument. 
Exemple : 

fichier mod_cte_math.f90 


module mod_cte_math 
implicit none 
real, parameter :: pi=3.141593 


end module mod _cte math 


fichier azre.f90 


program aire 
use mod_cte_math 
implicit none 
real :: r,s 


r=1.0 
S=pix*xr**2 


end program aire 


— permet aussi de passer des procédures en argu- 


ments, pour faire des fonctions de fonctions (voir 


section 5) 


— conseil : utiliser les modules autant que possible. 


4.5 Compilation séparée 


comment compiler un programme quand on a des sous- 


programmes externes situés dans des fichiers différents ? 
— compilation globale 


— utiliser la commande vue sec 2.9 


gfortran -o nom_executable liste_fichiers_fortrans 


— mais les modules doivent être compilés avant 
les sous-programmes utilisateurs. Ils doivent donc 


étre placés avant dans la liste 


— inconvénient : une modification d’un fichier en- 


traîne la re-compilation de tous les fichiers 


— compilation séparée : 2 phases successives 


1. compilation (option -c) : pour tout fichier fortran 
(ici fich.f90) faire 


gfortran -c fich.f90 


— création du fichier objet fich.o (de méme nom 


résultat : 


que le fichier fortran mais d’extension .0) 


— si fich.f90 contient le module fich, création du 
fichier fich.mod (de même nom que le module 
avec l’extension .mod). Contient Vinterface du 
module, nécessaire pour compiler les programmes 


utilisateurs de fich 


— attention : compiler les modules avant les sous- 


programmes utilisateurs 


2. édition de liens (option -o) : les fichiers objets (.0) 


sont liés ensemble pour faire l’exécutable 


gfortran -o nom_executable liste_fichiers_objets 


— intérêt : 
— on ne recompile que les fichiers modifiés et les fi- 
chiers dépendants, et l’on refait l’édition de liens 


— outil complémentaire efficace pour la compilation 


séparée de gros programmes (plusieurs dizaines de 


fichiers) : la commande make (voir la fin de cette 


section) 


— Exemple : programme fortran90 composé d’un pro- 
gramme principal prog qui utilise deux modules mod1 
et mod2. Le module mod1 utilise le module mod2. 


pr og exécutable 
) 4 
mas ® 


(3) | 


édition de liens 


mod1.0 
, , mod1.mod 


compilation (2) 


fichiers 
mod2.o iets 
mod2.mod 


(1) | 
prog.f90 mod2.f90 


=> la compilation séparée est : 
(1) compilation : 
(2) compilation : 
190 -c mod1.f90 
(3) compilation : 
f90 -c prog.f90 
(4) édition de liens : 


f90 -o prog prog.o mod1.o mod2.0 


= en cas de modification de prog, alors seules la 
compilation de prog.f90 puis l’édition de liens sont 


nécessaires (phases (3) et (4)) 


— utilisation de la commande : le fichier suivant 
(nommé “Makefile”) contient les dépendances entre 


fichiers et les actions à exécuter pour mettre à jour 


Vexécutable : 


# edition de liens 


prog : prog.o modi.o mod2.o 


<TAB> gfortran -o 


# compilation 

prog.o : prog.f90 
<TAB> gfortran -c 
modi.o : mod1.f90 
<TAB> gfortran -c 
mod2.o : mod2.f90 


<TAB> gfortran -c 


prog prog.o modi.o mod2.o0 


modi.o mod2.0o 
prog.f90 
mod2.0o 


modi.f90 


mod2.f90 


La commande unix utilise ce fichier pour ne 
a été modifié : il suffit de taper 


compiler que ce qui 


cette commande dans le dossier contenant le pro- 


gramme. 


Attention : la tabulation (visualisée ci-dessus par 
<TAB>) est indispensable. 


— on peut paramétrer ce “Makefile” pour utiliser des 


options de compilation différentes ou pour “nettoyer” 


le dossier : 


paramètres 
DEBUG= 
COMPILE=gfortran $ (DEBUG) 


# edition de liens 
prog : prog.o modi.o mod2.o 
<TAB> $(COMPILE) -o prog prog.o modi.o mod2.0 


<TAB> @echo compilation terminée 


# compilation 

prog.o : prog.f90 modi.o mod2.0 
<TAB> $(COMPILE) -c prog.f90 
modi.o : modi.f90 mod2.0o 

<TAB> $(COMPILE) -c mod1.f90 
mod2.0 : mod2.f90 

<TAB> $(COMPILE) -c mod2.f90 


# destruction des fichiers objets et modules 
detruire 


<TAB> rm -f *.0 *.mod 


— différents appels de : 


— | make |: compilation normale 


— | make OPT=-g | : compilation en mode debug 


— | make OPT=-O2]|: compilation optimisé de niveau 2 


— | make detruire | : efface les .o et les .mod 


— Voir |4, 1] pour des détails. 


— Remarque : pour les petits programmes, le make- 


file peut etre généré a la main, mais pour les pro- 


grammes plus longs, il peut étre utile de le générer 


automatiquement, par exemple avec le script perl 


mkmf|, voir : 


http://www.gfdl.noaa.gov/"vb/mkmf.html 


5 Utilisation avancée 


— les concepts présentés ici sont plus complexes que les 
précédents 
— a part le concept de précision des réels, les autres 


ne sont pas forcément nécessaires pour réaliser des 


programmes de calcul scientifique, sauf peut-être les 


pointeurs qui peuvent étre indispensables dans cer- 


tains cas 


— cependant, ces concepts peuvent permettre de sim- 


plifier l’écriture des programmes 


5.1 Types dérivés 


— structures de données renfermant plusieurs types 
différents 
— exemple : type personne contenant les nom, prénom, 


et age d’une personne : 


type personne 
character(len=30) :: nom,prenom 
integer :: age 

end type personne 


on peut alors déclarer des variables de ce type, comme 


pour n’importe quel type 


les différentes variables d’un type dérivé sont ap- 


pelées champs. On y accède avec le caractère %, et 
on peut utiliser le constructeur de même nom que le 
type pour les initialiser : 


etudianti=personne(’Cesar’,’ Jules’ ,57) 
etudiant2—etudiantil 
print x,etudiant2 


do i=1,2 
print,’ nom, prenom, age’ 
read x ,td(i)%nom ,td(i)%prenom ,td(i)%age 


end do 


print x ,td(1) 


l'affectation etudiant2=etudiant1 recopie tous les 
champs. 


5.2 Précision des réels 


— par défaut, le type real est dit ”simple précision”, 
c.-à-d. codé sur 4 octets (—32 bits), avec 6 chiffres 
significatifs : ceci est souvent insuffisant pour faire 
des calcul scientifique précis ; 

— il faut alors utiliser la représentation “double précision”, 
c.-à-d. codée sur 8 octets (—64 bits), avec 15 chiffres 
significatifs = 3 possibilités : 

— type double precision : 


double precision :: 


— type real avec paramètre kind : 


integer, parameter :: pr=kind(1.d0) 
real(pr) :: x 


(type réel de même précision que 1.40, c.-à-d. 
double précision) 
— ou 
integer , parameter :: pr=selected_real_kind (15,3) 


real(pr) :: x 


(type réel avec 15 chiffres significatifs et exposant 
a 3 chiffres, c.-a-d. double précision) 


— la 3°possibilité est la plus “portable” (donne la méme 
précision quelle que soit la machine) 


— les valeurs numériques s’écrivent différemment selon 


la précision choisie 


1 L lpr 


A, 32 4.32_pr 
4.32 x 107 * 4.32e—4_pr 


— attention à la cohérence des types dans les affecta- 


tions ou opérations : 


integer :: n 
double precision 
n 
.00000000000000 


{Í 
1.20000000000000 
1 
1 


. 20000004768372 


— par prudence : ne pas mélanger les entiers et réels 


— complexes : comme pour les réels : 


complex + À lee simple p. 
double complex :: z !--- double p. 
pr=kind (1.40) I! -- ou 
pr=selected_real_kind (15,3) 

complex (pr) o Z 1! --- double p. 


5.3 Fonctions de fonctions 


il s’agit simplement de créer des procédures dont les 


arguments peuvent être d’autres procédures 


exemple simple : écrire une fonction integrale qui 


pour toute fonction f et tout couple de réels a,b 


calcule une approximation de [° Ta) de 


pour cela : le fonction integrale doit être mise dans 
un module, et le type de la fonction f et de ses argu- 


ments doivent être déclarés dans un bloc interface 


à l’intérieur de la fonction integrale 


attention : dans le programme utilisant la fonction 


integrale, la fonction f passée en argument est : 

— soit une fonction intrinsèque (comme cos, log, 
exp, sqrt ...), il faut alors la déclarer avec l'at- 
tribut intrinsic ; 

— soit définie dans un module toto, il suffit alors 
de déclarer l’utilisation de ce module (use toto) ; 


— soit définie dans un sous-programme externe non 
placé dans un module, il faut alors la déclarer 
avec l’attribut external. 


— exemple : 


module mod_integrale 
implicit none 


contains 


real function integrale(f,a,b) 
implicit none 


l--- arguments 
real, intent(in) :: a,b 


!--- argument f : bloc interface 
interface 
real function f(x) 
real, intent(in) 
end function f 
end interface 


l--- yar locales 
integer :: n,i 
real :: h 


n=100 ; h=(b—a)/n 

integrale—0. 

do i—0,n-1 
integrale—integrale+hxf(ixh) 

end do 


end function integrale 


end module mod_integrale 


— utilisation de cette fonction dans un programme : 


program prog_integrale 


use mod_integrale 
use mod_creneau 
implicit none 


real :: yl,y2 
real , intrinsic :: exp 


yi=integrale(exp ,0.,1.) 
y2=integrale(creneau , —1.,1.) 
print*x,y1,y2 


end program prog_integrale 


module mod_creneau 
implicit none 
contains 
real function creneau(x) 
implicit none 
real, intent(in) :: x 


if (x<0.0.and.x>1.0) then 
creneau—(|. 
else 
creneau=l. 
end if 
end function creneau 
end module mod_creneau 


5.4 Interface générique 


— idée : on veut appliquer une même procédure à des 
arguments de type différents. Exemple : abs(x) ren- 
voie |x| si x est réel et Va? +b? si x = a + ib est 


complexe. 


principe : dans un module, on écrit plusieurs procédures 
adaptées à chaque type d’arguments, et on les re- 


groupe sous un même nom “générique”. C’est le com- 


pilateur qui décide quelle procédure utiliser selon le 
type des arguments. 


exemple : on veut une fonction exponentielle qui 


renvoie exp(x) si x est réel, et exp(z) = exp(a)(cos b+ 
isinb) si z = a + ib est complexe 


module mod_exponentielle 
implicit none 
interface exponentielle 
module procedure rexp, zexp 
end interface 
contains 
function rexp(x) result (y) 
implicit none 
real, intent(in) :: x 
real :: y 


y=exp (x) 
end function rexp 


function zexp(z) result(y) 
implicit none 
complex, intent(in) :: z 
complex :: y 
real :: a,b 


a=real(z) ; b=aimag(z) 


y=cmplx(cos(b),sin(b))*xexp(a) 
end function zexp 


end module mod_exponentielle 
— méthode : deux fonctions rexp et zexp stockées dans 


le module mod_exponentielle. La fonction générique 


exponentielle est créée avec le bloc interface et 


Vinstruction module procedure située au tout début 


du module. 


— ensuite, dans un programme principal, on utilise la 
même fonction exponentielle, que l’argument soit 


réel ou complexe : 


program generique 
use mod_exponentielle 
implicit none 


: 7—CHplx (0-3-1416) 


’ ,exponentielle(x) 


printx, ’exp(x) 
print x, ,exponentielle(z) 


end program generique 


5.5 Création de nouveaux opérateurs 


— exemple : on a une fonction rptens qui avec deux 


vecteurs réels vı et v en entrée renvoie le produit 


tensoriel w = v1 © ve : 


w=rptens(vi,v2) 


— on voudrait utiliser cette fonction sous forme d’un 


opérateur ptens : 


— on va utiliser un module, un bloc interface operator 


et l’instruction module procedure : 


module mod_ptens 


implicit none 
interface operator(.ptens.) 
module procedure rptens 
end interface 
contains 


function rptens(vi,v2) result(w) 
implicit none 
real, dimension(:), intent(in) :: vi,v2 
real, dimension(size(vi),size(v2)) :: w 
integer :: i,j 


do i=1,size(vi) 
do j=l,size(v1) 
w(i, j)=v1 (4) *v2 (5) 
end do 
end do 


end function rptens 


end module mod_ptens 


— les deux points autour de .ptens. sont obligatoires 


pour définir un opérateur 
— dans le programme principal utilisant le module 


mod_ptens, on pourra utiliser l’opérateur ptens ainsi: 
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— on peut associer plusieurs fonctions à un même opé- 
rateur (comme avec les interfaces génériques, sec- 
tion 5.4) : pour que .ptens. marche aussi avec les 
complexes, rajouter une fonction zptens dans le mo- 
dule, et ajouter le nom zptens à la suite de rptens 


dans l’instruction module procedure 


on peut aussi vouloir redéfinir un opérateur existant 


comme *,+,—,/ etc. Cela est possible, mais un peu 


plus compliqué (il faut que le nouvel opérateur et 


l’ancien diffèrent par le type de leurs arguments). 


Dans l’exemple précédent, on ne pourrait pas redéfinir 
x par l’opérateur ptens car * est déjà défini pour les 
tableaux. 
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— exemple de surcharge de + : pour concaténer des 


chaînes comme l’opérateur //, en supprimant des 


blancs. Remarquer que + garde sa signification pour 


les autres types. 


program surcharge 


use mod_add_chaine 
implicit none 


character(len—20) :: c1,c2,c 
él" bon” = C2” jour” 
e=ci7-c2 

print x,c 

printx,’hello da’+’ rling’ 
print *,5+3 


end program surcharge 


module mod_add_chaine 


implicit none 
interface operator(+) 
module procedure add_chaine 


end interface 
contains 
function add_chaine(c1,c2) result(c) 


implicit none 


character (len=>*), intent(in) :: c1,c2 
character (len=len(ci)+len(c2)) :: c 


c=trim(adjustl(ci))//trim(adjustl(c2)) 


end function add_chaine 


end module mod_add_chaine 


5.6 Ressources privées, publiques, semi- 


privées 


— quand un programme utilise un module, il a accès 
a toutes les variables du module, et tous ses sous- 
programmes. Cela peut parfois provoquer des conflits 
avec les propres variables du programme si les noms 


sont identiques 


exemple : un programmeur veut utiliser un module 
servant a calculer des normes de matrices avec la 
fonction norme pour son programme. Ce module con- 
tient beaucoup d’autres fonctions auxiliaires utilisées 
par la fonction norme. Le programmeur ne connait 
pas ces autres fonctions (il n’a accès qu’au fichier 
objet par exemple). Dans ce cas, les fonctions auxi- 
liaires ne devraient pas être visibles à l’extérieur du 
module, car sinon, l’utilisateur n’aurait par exemple 


pas le droit de créer des fonctions portant le même 
nom, alors qu’il ignore leur existence! 


on peut alors utiliser les fonctions private (privé), 
public et les attributs correspondants private et public 
pour limiter l’accès à certaines variables depuis l’ex- 


térieur du module 


— voir la référence [3] pour plus de détails. 


5.7 Pointeurs 


— la structure de pointeur est fondamentale en langage 
C et C++. En fortran 90, une structure similaire 


existe, mais n’est pas d’un intérêt évident pour le 


calcul scientifique 


dans la norme fortran 95, elle était indispensable 
pour définir des tableaux allouables à l’intérieur de 


types dérivés (comme par exemple pour définir des 


tableaux de tableaux). Mais depuis la norme fortran 


2000, cela n’est plus nécessaire. 


la notion de pointeur en fortran 90 ne sera donc pas 


étudiée plus en détail dans ce cours 
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Liste de quelques fonctions intrinsèques usuelles 
en Fortran 90 


adapté des cours de VIDRIS, disponibles sur 
http ://www.idris.Fr/data/cours/lang/f90/F90_cours.html 


On donne la signification de la commande, puis, dans un exemple, la syntaxe à utiliser, suivie 
après le signe + ’ de la valeur retournée en sortie. 


1 Fonctions intrinsèques sur les entiers, réels, complexes, 
extensibles aux tableaux 


— ABS : retourne la valeur absolue de son argument. Pour un complexe, retourne sa norme : 
1/22 + y?. 
ABS(-1) + 1 
ABS(—1.5) + 1.5; ABS((3.,4.)) + 5.0 


— ACOS : retourne l’arc cosinus en radians de son argument réel. 
ACOS(0.54030231) + 1.0 

— AIMAG : retourne la partie imaginaire du complexe passé en argument. 
AIMAG((2.,3.)) => 3.0 

— AINT : tronque le réel passé en argument 
AINT(2.783) +» 2.0; AINT(—2.783) => —2.0 

— ANINT : retourne, sous forme d’un réel, l’entier le plus proche du réel transmis. 


ANINT (2.783) — 3.0 
ANINT(—2.783) ++ —3.0 


— ASIN : retourne l’arc sinus en radians de son argument réel. 
ASIN(0.84147098) + 1.0 

— ATAN : retourne l’arc tangente en radians de son argument réel. 
ATAN(1.5574077) + 1.0 

— CEILING : retourne l’entier immédiatement supérieur au réel transmis en argument. 


CEILING(3.7) + 4 
CEILING(—3.7) + -3 


— CMPLX : retourne un complexe dont les parties réelle et imaginaire sont transmises en argument. 
CMPLX(—3.) 1% —3.0+0.i; CMPLX(2,4.) + 2.0+4.0i 

— CONJG : retourne le complexe conjugué de celui passé en argument. 
CONJG((—3.0,4.0)) => —3.0—-4.0i 

— COS : retourne le cosinus de l’angle passé en argument (exprimé en radians). 


COS(1.0) + 0.54030231 
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COSH : retourne le cosinus hyperbolique. 
COSH(1.0) + 1.5430806 


DBLE : convertit en double précision l’argument transmis. 
EXP : retourne l’exponentiel de l’argument transmis. 


EXP(1.0) + 2.7182818 

FLOOR : retourne l’entier immédiatement inférieur au réel transmis en argument. 
FLOOR(3.7) — 3 ; FLOOR(—3.7) — —4 

INT : convertit en entier l’argument transmis. 


INT (327) <3 
INT (9.1/4.0) + 2 


LOG : retourne le logarithme népérien de l’argument transmis. 


LOG(2.7182818) — 1.0 
LOG(10.0) + 2.3025851 


LOG10 : retourne le logarithme décimal de l’argument transmis. 


LOGO CIO 0j = 210 
LOG10(10.E10) — 11.0 


MAX : retourne le maximum des nombres passés en argument. 
MAX(—9.0,7.0,2.0) + 7.0 

MIN : retourne le minimum des nombres passés en argument. 
MIN(—9.0,7.0,2.0) + —9.0 

MOD : retourne le reste de la division effectuée à l’aide des deux arguments fournis. 


MOD(3.0,2.0) ++ 1.0 
MOD(—8,5) > —3 


REAL : convertit en réel l’argument transmis. 
REAL(3) + 3.0 


SIGN : retourne le nombre dont la valeur absolue est celle du premier argument et le signe celui 
du deuxième. 


SIGN(—3.0,2.0) + 3.0 

SIN : retourne le sinus de l’angle passé en argument (exprimé en radians). 
SIN(1.0) + 0.84147098 

SINH : retourne le sinus hyperbolique. 
SINH(1.0) + 1.1752012 

SQRT : retourne la racine carré de son argument. 
SQRT(5.0) + 2.236068010 

TAN : retourne la tangente de l’angle passé en argument (exprimé en radians). 
TAN(1.0) + 1.5574077 

TANH : retourne la tangente hyperbolique. 
TANH(1.0) + 0.76159416 


109 


2 Fonctions intrinsèques pour le temps de calcul 


— CPU_TIME(t) (Norme 95) retourne dans le réel t le temps CPU en secondes (ou réel < 0 
si indisponible). Par différence entre deux appels, il permet d'évaluer la consommation CPU 
d’une section de code. 


call cpu_time(ti) 


call cpu_time(t2) 
print,’ temps CPU_—’ ,t2-t1 


3 Fonctions intrinsèques pour les chaînes de caractères 
— ADJUSTL : cadre à gauche la chaîne passée en argument : supprime les blancs en tête ; complète 
à droite par des blancs. 
ADJUSTL(”’... Fortran’) + ‘Fortran? 


— ADJUSTR : cadre à droite la chaîne passée en argument : supprime les blancs en fin; complète 
à gauche par des blancs. 


ADIUSTR( Fortransass les is Fortran 
— LEN : retourne la longueur de la chaîne de caractères transmise en argument. 
CHARACTER(len—10) CH; LEN(CH) ++ 10 


— LEN_TRIM: retourne la longueur de la chaîne de caractères transmise en argument sans considérer 
les blancs de fin. 


LEN_TRIM( „FORTRAN u2’) > 9 
LEN_TRIM(’uuu’) > 0 


— LGE : compare les deux chaînes de caractères transmises en argument : retourne .true. si la 
première chaîne est supérieure ou égale à la deuxième, .false. sinon. 


LGE ( ’ MANET’ , MONET’) + false. 
LGE ( MANET Edouard’ , MANET’) + .true. 


— LGT : compare les deux chaînes de caractères transmises en argument : retourne .true. si la 
première chaîne est supérieure strictement à la deuxième, .false. sinon. 


LGT ( MANET’, MANET’) ++ .false. 


— LLE : compare les deux chaînes de caractères transmises en argument : retourne .true. si la 
première chaîne est inférieure ou égale à la deuxième, .false. sinon. 


LLE ( ’MANET’ , MONET’) + .true. 
LLE( ’MANET’ , MANET’) + true. 


— LLT : compare les deux chaînes de caractères transmises en argument : retourne .true. si la 
première chaîne est inférieure strictement à la deuxième, .false. sinon. 


LLT( MANET”’ , MANET’) + false. 

— TRIM : retourne la chaîne de caractères transmise débarrassée de ses blancs de fin. 
TRIM( PICASSO... 7) + ° PICASSO’ 

— // : opérateur pour concaténer deux chaînes. 


‘bon’ // ‘jour’ + ‘bonjour’ 
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4 Fonctions et subroutines intrinsèques tableaux 


4.1 interrogation 


— SIZE(array[,dim]) retourne la taille (ou l’étendue de la dimension indiquée via dim) du 
tableau passé en argument. 
integer, dimension(—2:27,0:49) :: t 
SIZE(t) + 1500 
SIZE(t,dim=1) — 30 
— MAXLOC(tab) et MINLOC(tab) retournent, pour le tableau tab, l'emplacement de l’élément 
maximum/minimum (si le tableau est multidimensionel, le résultat est un vecteur dont les 
éléments sont les indices de l’élément dans tab) 
integer, dimension(—2:27,0:49) :: t 
integer, dimension(1:2) :: ind_elem_max, ind_elem_min 
ind_elem_max=maxloc(t) 
ind_elem_min=minloc(t) 


4.2 Réduction 


Selon que DIM est absent ou présent, toutes ces fonctions retournent soit un scalaire soit un 
tableau de rang n-1 en désignant par n le rang du tableau passé en premier argument. 
La plupart de ces expressions acceptent un “masque” en paramètre, qui permet de n’appliquer 
les fonctions qu’à certains éléments (voir le cours de l’IDRIS de Fortran 95) 
— COUNT(verif) : verif est une expression logique sur deux tableaux conformes, count compte 
le nombre de fois où cette expression est vérifiée par deux éléments 
SE: 19, Be" 
Soit A | | 


Décompte global des valeurs de A > 2 
COUNT(A>3) = 4 


— MAXVAL(array) et MINVAL (array) 


ae ee 
Soit A | 

12 4 6 
maxval(A) — 6 
minval(A) + 1 


— PRODUCT (tab) renvoie le produit de tous les éléments du tableau tab. 
PRODUCT ((/ 2,5,-6 /)) + —60 
— SUM(tab) renvoie la somme de tous les éléments du tableau tab. 


SUM((/ 2,5,-6 /)) + 1 


4.3 multiplication 


— DOT_PRODUCT(vector_a,vector_b) retourne le produit scalaire des deux vecteurs passés en 
argument (produit hermitien si vector_a est complexe) 


111 


vi =(/ 2,—3,—1 /) 3 v2 =(/ 6,3,3 /) 
DOT_PRODUCT(v1,v2) + 0 


— MATMUL (matrix_a,matrix_b) effectue le produit matriciel de deux matrices ou d’une matrice 
: l 3 —6 —1 2 
et d’un vecteur passés en argument. Soit la matrice A = ( 2 3, 1 ) et le vecteur v = (<) 


MATMUL(A,V) 


4.4 transformation 


— TRANSPOSE(mat) renvoie la transposée de la matrice mat (tableau bidimensionnel) 
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